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We review our results for the QCD phase diagram at baryonic chemical potential 
u B < wT. Our simulations are performed with an imaginary chemical potential ui 
for which the fermion determinant is positive. For 2 flavors of staggered quarks, we 
map out the phase diagram and identify the pseudo-critical temperature T c (fij). 
' For fJ-i/T < 7r/3, this is an analytic function, whose Taylor expansion is found 

to converge rapidly, with truncation errors far smaller than statistical ones. The 
truncated series may then be continued to real fi, yielding the corresponding phase 
diagram for /is<500 MeV. This approach provides control over systematics and 
|^Nj ' avoids reweighting. We outline our strategy to find the (2+l)-flavor critical point. 

1. Introduction 

Substantial progress has been accomplished over the last two years toward 
the numerical determination of QCD properties at finite density. This 
progress has not come through a solution of the "sign problem" , which 
prohibits standard Monte Carlo simulations because the Dirac determinant 
is no longer real positive at non-zero chemical potential [i. The standard 
approach to the sign problem, which requires statistics growing exponen- 
tially with the volume and the chemical potential, is still current. Rather, 
progress has come through a change to a more pragmatic attitude. It makes 
sense to explore the region of small chemical potential, with methods whose 
failure at large fi has become acceptable. This is because information so ob- 
tained is phcnomenologically important, especially for heavy-ion collisions 
where the net quark density remains small. 
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Three approaches are currently pursued, and reviewed in these proceed- 
ings: (i) a two-parameter reweighting, in T and fi, of \i = simulations 2 ; 
(ii) a Taylor expansion of (i) truncated to its lowest non- trivial order 3 ; (Hi) 
a study at imaginary /i, followed by analytic continuation of a truncated 
Taylor series to real ll. (i) is limited to small volumes because of the sign 
problem; (ii) is in addition restricted to smaller chemical potentials, and 
the systematic error due to the Taylor truncation is unknown. We have 
adopted (Hi): it can be used for arbitrarily large volumes since it has no 
sign problem, and systematic errors coming from analytic continuation can 
be controlled. We review this approach here, applied to the determination 
of the pseudo-critical line T c (p) between confinement and deconfinement 1 . 

2. QCD at real and imaginary chemical potential 

The grand canonical partition function Z(V,/i,T) = Tr (^q-(h~vQ)/t^ can 
be considered for a complex chemical potential \i — llr + i/ij. Two general 
properties suffice to constrain the phase structure as a function of 

- Z is an even function of ji: Z(jl) — Z(—p,), where fx = ll/T . 

- A non-periodic gauge transformation, which rotates the Polyakov loop by 
a center element but leaves Z unchanged, is equivalent to a shift in lii 4 : 

Z(fx R ,fn) = Z(fx R ,fx I + 2n/N). (1) 

For QCD (N = 3), these two properties lead to Z(3) transitions at critical 
values of the imaginary chemical potential, flj = ^ (n+ ^), separating 
regions of parameter space where the Polyakov loop angle (ip) falls in dif- 
ferent Z(3) sectors. Perturbative and strong coupling considerations led 
the authors of 4 to predict at pff a first-order phase transition at high tem- 
perature (deconfined phase) and a crossover at low temperature (confined 
phase). The resulting phase diagram in the (ni,T) plane is periodic, and 
depicted in Fig. 1. 

The non-analyticity closest to the origin, at jlj = ^, limits the prospects 
of analytic continuation. Our strategy, following 5 , consists of fitting by 
a Taylor series in jij observables measured at (fin — 0,fti ^ 0), then 
continuing the truncated Taylor series to real fi. Convergence of the Taylor 
expansion can only be checked for lii < |T, i.e. /zs<500 MeV. 

3. Analyticity of the (pseudo-) critical line 

The phase transition or crossover between confinement and deconfinement, 
as (3 (or T) is varied while keeping /i fixed, is characterized by a peak of 
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Figure 1. Schematic phase diagram in the (/nj,T) plane. The solid line marks a first- 
order Z3 transition. 



the susceptibility x = VN t ((O — (O)) 2 ) for observable O. The critical line 
/3 c (/x) is thus defined implicitly via 



dp 



d 2 X 

= — — 

dp 2 



< 0. 



(2) 



In a finite volume V, xiPiN i s analytic, and from the implicit function 
theorem j3 c (fj,) also is. Moreover, it is an even function of \i just like Xi so 
that it can be expanded in powers of fi 2 : 



p c { N =/3 c (/i = 0) + ^c„(a M ) 2 ". 



(3) 



71=1 



Analytic continuation between real and imaginary \i (for < §T) is 
accomplished by flipping the sign of [i 2 . 



4. Numerical results for two light flavors 

We have studied QCD with 2 flavors of staggered fermions (8 3 x 4 lattice, 
quark mass am q — 0.025, R-algorithm) . The lattice spacing is a ~ 0.3 fm, 
and sizeable corrections can be expected in the continuum extrapolation. 

The predicted nature of the Z(3) transition at (a/i/) c = 7r/12 is 
confirmed, and the phase diagram is that of Fig. 1. The confinement- 
deconfinement transition line is obtained from the peaks of the plaquette, 
quark condensate and Polyakov loop susceptibilities, which give consistent 
values P c (ani) (see Fig. 2). These values are then fitted by a Taylor series 
in {afij) 2 over the interval am € [0, n/12]. The data are well described by a 
linear fit, and agreement hardly improves with a quadratic fit (see Fig. 3a). 



1, 2008 20:19 WSPC/Trim Size: 9in x 6in for Proceedings 



myproc34 




The quadratic coefficient oc (a/i/) 4 is zero within errors. Therefore, we can 
safely approximate the critical line as 

/3 c (a/ir) = 5.2865(18) + 0.596(40)(a^) 2 . (4) 

We translate into physical units using the perturbative two-loop /3-function, 
which suffices for our present accuracy. This yields, as a function of the 
baryonic chemical potential /is 

5^L = l- 0.00563(38) (f) " , (5, 

using T c {fi = 0) = 173(8) MeV, while the next-order term 0((/i S /T) 4 ) is 
statistically insignificant up to /xg ~ 500 MeV (see Fig. 3b). Similar results 
have been obtained for 4 flavors 6 . 



5. How to find the (2+l)-flavor critical endpoint 

In the case of 3 flavors (degenerate or not), the deconfinement "transition" 
changes from crossover to first-order at a critical point (T c ,/x c ), which can 
be identified by a study of cumulant ratios or Lee- Yang zeroes. As the 
quark mass m q is varied, the critical point describes an analytic curve 
T c (/i) (Fig. 4), which can be Taylor expanded. On the lattice: 

/3 C ( M ) = = 0) + £>nM 2n ■ (6) 

n=l 

As in Section 4, data can be obtained at imaginary /i, then analytically 
continued to real fx by flipping the sign of /i 2 . 
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Figure 3. (left) Phase diagram in the ((a/j,j) 2 , j3) plane. The rightmost point lies beyond 
the Zz transition, and is therefore excluded from the fit. (right) Analytically continued 
phase diagram in the (fig,T) plane. 
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Figure 4. Critical lines in the (T, fi) plane for different quark masses m q . The bold 
upward parabolic curve characterizes second-order transition points, separating the 
crossover and the first-order regimes. It can be analytically continued to negative /x 2 . 
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